Feshbach-resonant Raman photoassociation in a Bose-Einstein condensate 
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We model the formation of stable heteronuclear molecules via pulsed Raman photoassociation of 
a two-component Bose-Einstein condensate near a strong Feshbach resonance, for both counterintu- 
itive and intuitive pulse sequencing. Compared to lasers alone, weak Raman photoassociation is en- 
hanced by as much as a factor of ten (five) for a counterintuitive (intuitive) pulse sequence, whereas 
strong Raman photoassociation is barely enhanced at all-regardless of pulse sequence. Stronger 
intra-atom, molecule, or atom-molecule collisions lead to an expected decrease in conversion effi- 
ciency, but stronger ambient inter-atom collisions lead to an unexpected increase in the efficiency of 
stable molecule production. Numerical results agree reasonably with an analytical approximation. 



PACS numbers: 03.75.Fi,05.30.Fk,32.80.Wr 

I. INTRODUCTION 

Not long after a Bose-Einstein condensate of atoms 
was created [T], the race began to reach the molecu- 
lar milestone, and thereby further enable fundamental 
studies [2J [3], practical applications g], and proxy in- 
vestigations [5J [5J . The problem is that the laser cooling 
techniques that enable atomic condensation are difficult- 
though not impossible-to apply to molecules [Jj, while 
Stark deceleration [5J and buffer gas cooling [5] have yet 
to reach quantum degeneracy. Alternatively, association 
of atoms into molecules-driven by a laser [10] or mag- 
netic field [H]-is highly efficient at quantum degenerate 
phase space densities [12] , and a fundamental coherence 
means that a condensate of atoms could be associated 
into a condensate of molecules [13]. Indeed, high effi- 
ciency P3] [T3] and atom-molecule coherence [TH1 E] have 
been demonstrated in both photoassociation [T3J [T7J and 
magnetoassociation |15L I16j . and quantum degeneracy 
has been achieved on short time scales with magnetoasso- 
ciation |18j . However, production of a long-lived conden- 
sate of molecules in the absolute ground state is ham- 
pered by collisions between particles and limited laser 
intensity [19] , and pre-loading an atomic condensate into 
an optical lattice [20] was therefore proposed to mitigate 
collisional effects in association [21 . So far, a quantum 
degenerate gas of molecule in the absolute ground state 
has been produced by several groups [22] , but melting of 
the lattice to create a bulk condensate of stable molecules 
has yet to be achieved. 

Here we explore a different route to forming a bulk 
condensate of molecules in the absolute ground state. In 
particular, combining the photoassociation and Feshbach 
resonances has been shown to enhance the photoassoci- 
ation rate constant [231423] of a Bose-Einstein conden- 
sate by an order of magnitude [2H [25] , due to construc- 
tive quantum interference between molecules formed by 
direct photoassociation and molecules formed by pho- 
toassociation via the Feshbach molecular state [25] . We 



therefore consider bound-free-bound-ground transitions 
[Fig. [lja)], whereby a magnetic field converts atom 
pairs into vibrationally-excited Feshbach molecules, a 
photoassociation pump laser converts atoms pairs into 
electronically-excited molecules, and a secondary dump 
laser converts the photoassociated molecules into stable 
molecules [2SJ [27] . The lasers are pulsed to best avoid 
irreversible losses, and we consider counterintuitive se- 
quencing where the dump pulse precedes the pump pulse, 
as well as intuitive sequencing where the pump pulse pre- 
cedes the dump pulse. The model includes relevant elas- 
tic collisions between particles, dissociative decay of the 
Feshbach molecules, as well as spontaneous and dissocia- 
tive decay of the electronically-excited photoassociation 
molecules. A quasicontinuum model explicitly includes 
dissociation and Feshbach molecules, whereas a resonant- 
interaction model treats the dissociation continuum and 
Feshbach molecules as virtual-leading to a magnetically- 
tunable photoassociation coupling and dissociation rate, 
in addition to the usual magnetically-tunable collisional 
interaction. Both the full and resonant-interaction mod- 
els contrast with a previous model [25] where the Fesh- 
bach resonance was accounted for merely with a magnet- 
ically tunable collisional interaction. 

The question is whether Feshbach enhancement of pri- 
mary photoassociation into excited molecules will carry 
over into enhancement of Raman photoassociation into 
stable molecules. On one hand hand, conversion in the 
counterintuitive scheme is presently hobbled by a com- 
bination of collisions and laser intensity limited by con- 
densate size [T!5] , and Feshbach-enhancement of the pho- 
toassociation coupling delivers strong coupling at low 
intensity which, in turn, could mitigate collisions. On 
the other hand, only an odd number of intermediate 
levels can form a dark state |29| . and a counterintu- 
itive scheme involving Feshbach and photoassociation 
molecules should therefore fall short of its trademark unit 
efficiency. Then again, it is possible that the dissoci- 
ation continuum could act like a distinct intermediate 
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FIG. 1: Few-level diagram of Raman photoassociation near 
a Feshbach resonance. In a full model (a), the Feshbach and 
photoassociation molecules share the same quasicontinuum 
(QC) of dissociation states. In a four-level approximation 
(b), the virtual quasicontinuum shifts the photoassociation 
and Feshbach detunings, and effectively couples transitions 
between the two excited molecular states. In the three-level 
approximation (c), the virtual Feshbach molecular state shifts 
the photoassociation detuning, and effectively renormalizes 
the photoassociation coupling and detuning. 



level, as it has on a separate occasion [3U], enabling a 
fully effective counterintuitive scheme. Finally, if a coun- 
terintuitive scheme is less than perfect, the question of 
enhancement over laser alone still remains and, more- 
over, whether or not it outperforms the intuitive scheme. 
All told, whereas strong photoassociation is already sat- 
urated and thus essentially un-enhanced, the Feshbach 
resonance enhances weak Raman photoassociation in a 
condensate and, somewhat surprisingly, the enhancement 
is stronger for stronger ambient inter-atomic collisions. 



Our work is outlined as follows. In Sec.[TTJ we introduce 
the full model that includes elastic collisions between par- 
ticles, a shared dissociation quasicontinuum for the Fes- 
hbach and photoassociation molecules, as well as sponta- 
neous decay of the photoassociation molecules. Also, we 
develop the resonant-interaction model based on virtual 
continuum and Feshbach states, wherein the photoasso- 
ciation interaction, the photodissociation rate, and the 
s-wave collisional interaction between atoms are magnet- 
ically tunable. Section |III| briefly reviews the computa- 
tional algorithm and provides parameters for numerical 
experiments. Section |TV| reports results for the efficiency 
of stable molecule production as a function of magnetic 
field, and compares these results against the two-photon 
rate constant. Finally, a summary is given in Section [Vj 



II. QUASICONTINUUM AND 
RESONANT-INTERACTION MODELS 

We focus on a two-component condensate because the 
resulting heteronuclear molecules are of significant re- 
cent interest for their dipolar properties [2H1| (see also 
Refs. [3T]). Nevertheless, any dipolar interaction is con- 
sidered relevant only after the molecules are formed, and 
is not accounted for in the association process. We also 
expect the results to apply to homonuclcar systems. 

In the few- level description illustrated in Fig.[lja), we 
consider N% (N2) atoms of species 1 (2) that have Bose 
condensed into the state |0i) (IO2)), say, the plane-wave 
state with zero momentum Hk = 0, which are repre- 
sented in Fig. [lja) as a single state |0) = 1 0i ) | O2 ) ■ A 
magnetic field tuned nearby a Feshbach resonance then 
couples two atoms in |0), one from each condensate, to a 
vibrationally-excited molecule in the state |2). Addition- 
ally, a photoassociation pump laser couples the same two 
atoms to an electronically-excited molecule in the state 
1 3), and a secondary dump laser couples the molecule in 
1 3) to a molecule in the absolute ground state |4). In 
the full model of Fig. [lja), the Feshbach and photoas- 
sociation molecules dissociate into non-condensate atom 
pairs that occupy one of a quasi-continuum of states, say, 
plane- wave states of momentum ±/ik 7^ 0. 

In second-quantized notation, the Hamiltonian corre- 
sponding to Fig. [lja) is 

H = H M +Hp + H D + H c , (1) 
where the contribution due to the magnetic field is 

— - = loqc'c + ^ «k(c t ak,ia-k,2 + a_k,2 a k,i c )» ( 2 ) 

k 

the contribution due to the photoassociation pump laser 
is 

— = S tfb + ^ ^k(& t ak,ia-k,2 + a -k,2 a k,i & )' ( 3 ) 

k 
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the contribution due to the secondary dump laser is 



equations of motion are given by 



(4) 



and the contribution due to s-wave collisions in the con- 
densates is 



He 
h 
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Here atoms in the ith condensate are represented by 
a i = a i,0: atoms with momentum Kk by a^, Fesh- 
bach molecular condensate by Co = c, photoassociation 
molecules by bo — b, and stable molecular condensate by 
go = g. The detuning of the magnetic field from the Fesh- 
bach resonance is oj , spontaneous decay of the Feshbach 
molecules [32] is neglected, the one-photon laser detuning 
is 8 Q = 9t[5o]) the spontaneous decay rate for the photoas- 
sociation molecule is T s — 23[<5o], and the two-photon 
detuning is A. The magnetic-field coupling between the 
atoms and the Feshbach molecules is «k = Qf/M,k, the 
pump-laser coupling between the atoms and the pho- 
toassociation molecules is 57k — ^/p,ki and the dump- 
laser coupling between the photoassociation and stable 
molecules is \- The momentum dependence of the Fes- 
hbach and photoassociation couplings are contained in 
/m,1c and /p.ki respectively, where /k=o = 0- Finally, the 
strength of collisions is determined by Ay, which is de- 
termined by the s-wave scattering length. Compared to 
the spontaneous decay rate, elastic collision involving pri- 
mary photoassociation molecules are neglected, and we 
also neglect vibrational relaxation [33J of the photoasso- 
ciation and Feshbach molecules. 

The quasicontinuum mean-field model is derived from 
a c-number approximation to the Heisenberg equations, 
ihx — [x,H], with x is the relevant operator, which 
generally works best for N > 100 [34] . Dissociation 
of Feshbach and photoassociation molecules into non- 
condensate atoms pairs is accounted for with the opera- 
tor ak,ia_k,2, and the corresponding c-number amplitude 
A^ = (ak,ifl-k,2)- Lastly, the quasicontinuum in momen- 
tum is converted into a continuum in frequency according 
to Ek -> A7(47r 2 wp /2 ) J de, where He = h 2 k 2 /(2p) is the 
kinetic energy and uj p = Hp 2 / 3 /(2p) is the characteristic 
frequency for a dissociated pair, with p the total particle 
density and p the reduced atomic mass. The resulting 
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Defining A^j — \pXij and Ai = j — pXij, the re- 
spective mean-field shifts are Ai = An|ai| 2 + Ai 2 |a 2 | 2 , 
A 2 = A 22 |a 2 | 2 + Ai 2 |ai| 2 , A c = A 00 |c| 2 + A 03 | 5 | 2 , and 
A 9 = A 33 | 52 | 2 + A 31 | ai | 2 + A 32 |a 2 | 2 + A 30 |c| 2 . Also, 
Cm = a/(4ir 2 ujp^ 2 ) is the magnetodissociation coupling 
and £p = f2/(47r 2 cjp^ 2 ) is the photodissociation coupling. 

The resonant-interaction model is derived by first 
treating the dissociated pair amplitude adiabatically 
(A = 0), which is equivalent to the limit of weakly-bound 
molecules 12 5|, and leads to the effective four- level sys- 
tem [Fig. [Tj(b)] with mean-field equations of motion 
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The virtual continuum leads to an effective coupling be- 
tween the Feshbach and photoassociation molecules [55] 
of strength 
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Additionally, there are real and imaginary shifts for each 
molecular detuning: uj = cj — cr m ~ *7m/2 and 6 = S — 
(Jp-i^p/2, where ct M (P) = ^Pm(p)] an d 7 l = 3[S M(P) ] 
with 
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(9a) 
(9b) 



The real shift is the well known result of coupling a bound 
state to a continuum [35] . and the imaginary shift is 
the dissociation rate. Of course, here the Feshbach shift 
um is static, while the photoassociation ap shift is tran- 
sient, and both are treated as implicit in the detuning 
rather than explicitly. Next, we arrive at the resonant- 
interaction model by treating the Feshbach amplitude 
adiabatically, which is the limit of large detuning of the 
magnetic field from the Feshbach resonance, and leads to 
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an effective three-level system, 

id\ — A\ai + QRCi^b (10a) 

id-x = A' 2 a2 + flRa\b (10b) 

ib = vb + Q^a\a2 + X9 (10°) 

ig = -(A-A g )g + X b. (lOd) 

As detailed previously [2"3"H2"5] . in the resonant- 
interaction model the Feshbach resonance effectively 
modifies the photoassociation interaction Qr = Q, — olk/oj 
and detuning v = vr — iTr/2, where the effective de- 
tuning is vr = S — k 2 /lu and the effective decay rate 
is T R = T + (k 2 /^ 2 )7 M with T = T s + j P . Depend- 
ing on the sign of the Feshbach detuning, the modified 
photoassociation coupling can be greater than the un- 
modified coupling, zero, or less than the unmodified cou- 
pling. Similarly, the real part of the resonant contribu- 
tion to the photoassociation detuning, k 2 /uj, can pro- 
duce a redshift, no shift, or a blueshift, as opposed to 
the ambient shift, <7p, which is strictly to the red |35j . 
Moreover, the resonant contribution to the decay rate, 
(k 2 /lu 2 )jm , leads to decay that depends on magnetic 
field and diverges near the Feshbach resonance. Fi- 
nally, we find the usual resonant collisional interaction, 
Ai = An|ai| 2 + A fl |a 2 | 2 and A 2 = A 22 |a 2 | 2 + A fl |ai| 2 , 
where A^ = Ai 2 — a 2 /2u>. In one-photon transitions [25] . 
collisions-resonant or otherwise-are neglected compared 
to the spontaneous decay rate, but here the timescale for 
conversion is long enough, especially in the counterintu- 
itive scheme, that collisions become relevant. 

III. NUMERICAL DETAILS 

For either the quasicontinuum model (|6| or resonant- 
interaction model ( |10p , the mean- field equations of mo- 
tion are written in matrix form, ihtp = Mip, where the 
nonlinear matrix M(ip) depends on the wavevector ip. 
Given the solution ip(t), the solution t(>(t + di) is then ap- 
proximated as i/)(t+dt) = M_M^ l tp(t), where M± = 1± 
i(dt/2)M, Given the solution to M + </>(t) = ip(t), the de- 
sired solution is simply ip(t + dt) = M_cf>(t). The nonlin- 
earity is accounted for with a two-step predictor-corrector 
method, wherein the prediction is tp p (t + dt) = M_<fi(t), 
and the correction is then ip(t+dt) — M_ </>(£) with 0(t) = 
M+rp(t), M± = M±(4>) and V^) = [ip p (t)+ip(t)]/2. This 
algorithm has the advantage that it can be made to scale 
linearly with the number of states [551 OH] , which is im- 
portant in the quasicontinuum model since we fix the 
number of quasicontinuum states to 10 6 . 

Considering the parameter values, the strength of the 
free-bound association couplings are measured relative to 
the characteristic frequency uj p [301 EH] and, for consis- 
tency, so are the collisional couplings. Hence, we con- 
sider a strong Feshbach coupling a 3> ui p , as well as both 
strong (Qq 3> w p) an d weak (f2 < oj p ) photoassociation, 
respectively. The laser pulses are modeled as Gaussian, 



n = n exp[-(i-7i) 2 /T 2 ] and x = Xo exp[-(t-T 2 ) 2 /r 2 ], 
so that k = hq cxp[— [(t — Ti) 2 /r 2 ]. Finally, collisions are 
pre-mitigated to a certain extent with low density |28j 
p = I0 12 cm~ 3 , and we assign elastic collisional couplings 
that range from weak to moderate. 

In particular, the Feshbach coupling is a — 134.8, 
the weak (strong) photoassociation pump coupling is 
fio = 15-4 (154), and the spontaneous decay rate is 
r s = 41S1 . Here strong photoassociation corresponds to 
a laser set at the saturation intensity and, since 
the weak coupling corresponds to an intensity two or- 
ders of magnitude below the saturation intensity |37l - 
HDj . In the counterintuitive scheme, the dump coupling 
is xo = 500o, but the intuitive scheme is more efficient 
for xo — S^o- The respective coupling to the dissoci- 
ation continuum is then determined by the Lorentzian 
fi = 1/(1 + e 2 //3 2 ), which is in turn determined by the 
molecular size /3, = H/mL 2 . For magnetoassociation we 
choose a point particle, Lm = P and for photoas- 
sociation we choose a typical size, Lp = I30ao, where 
ao is the Bohr radius. For the collisional couplings, 
we choose values for Ay such that A o = A 30 = 0.8, 
An = 5.1 x 10" 3 , A 22 = 2.1 x 10~ 2 , A12 = 4 x 10~ 2 , 
A 33 = 8.1 xl0~ 3 , A 3 i = 4.3xl0~ 2 , and A 32 = 1.9 xl0~ 2 . 
Finally, in the tunable collision model, the magnetodisso- 
ciation rate is Jm/^s — 7.4 x 10~ 2 , the photodissociation 
rate for strong (weak) coupling is [37 "fp/T s = 1 (1/100), 
and the peak cross-molecular coupling is k — 3.9 x 10 3 
(3.9 x 10 4 ) for weak (strong) photoassociation. 

To give the counterintuitive scheme the best chance for 
success, i.e., slow enough to foster adiabaticity but fast 
enough to outrun ambient collisions, we set the pulse 
width according to fio r = 5x 10 3 . Also, to ease the nu- 
merical overhead we only optimize the one-photon detun- 
ing S, but we fix the pulse delay to D = T\ — T 2 = — 2r, 
and we also fix the two-photon detuning to resonance 
(A = 0) in the weak case and Stark-shifted resonance 
(A = Oq/2<5) in the strong case. For the intuitive scheme, 
coincident pulses (D = 0) are generally optimal for all 
magnetic fields, both the pulse width and one-photon de- 
tuning are optimized at each magnetic fields, but the two- 
photon detuning is again fixed to resonance in the weak 
case and Stark-shifted resonance in the strong case. Al- 
though we do not go into details, for both pulse schemes 
the optimized one-photon detuning is consistent with the 
expected [231 dispersive-like behavior. 



IV. FESHBACH ENHANCEMENT 

This section reports the results of numerical experi- 
ments, summarized in Fig. [2] The node below resonance 
arises from destructive interference between direct pho- 
toassociation and indirect photoassociation occurring via 
the Feshbach molecular state [25 . A second node appears 
on resonance, and in the full model this node is due to 
the absence of an atom-molecule dark state that includes 
both excited molecules, while in the effective model is due 
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to divergence of the magnetically-tunable photodissoci- 
ation term. Regarding the magnitude of peak enhance- 
ment, for weak coupling the magnetic field enhances both 
schemes by about a factor of five, and the counterintu- 
itive scheme outperforms the intuitive scheme by about 
10% in the full model, and the two are roughly tied in 
the resonant-interaction model. Strong photoassociation 
is enhanced very little, to about about 80%, regardless of 
pulse order. The agreement with the full quasicontinuum 
model is reasonable given the simplicity of the effective 
model. It is worth noting that a model where the pri- 
mary photoassociation molecule is virtual instead of the 
Feshbach molecule improves the agreement between the 
two models near magnetic resonance, but worsens the 
off- resonant disagreement. 

Focusing on weak photoassociation, improvements to 
the counterintuitive scheme can be made by increas- 
ing the pulse area from fi T = 5 x 10 5 to 2.2 x 10 5 
(6 x 10 5 ) in the resonant-interaction (quasicontinuum) 
model, whereby the maximally-enhanced counterintu- 
itive efficiency improves from 30% to nearly 50% (76%). 
In the resonant-interaction model, further improvements 
can also be made by setting the pulse area according 
to flu instead of Q, but this would require a change in 
the dump coupling to \ — 50Qr, in order to satisfy the 
conditions for adiabatic following (specifically, x ^> 
for t — > — oo and fi 3> % for t — > oo), and therefore 
would require an impractical amount of dump laser in- 
tensity as D,r diverges on resonance. More importantly, 
this improvement is misleading since a larger peak dump 
pulse does not improve conversion in the full quasicontin- 
uum model, which we attribute to the absence of a dark 
state that includes both photoassociation and Feshbach 
molecules. 

The dip near resonance and the off-resonant peak to- 
gether indicate that Feshbach enhancement does not 
tell the whole story, since Feshbach enhancement peaks 
near magnetic resonance [55] • In the full model, 
cross-coupling between the photoassociation and Fesh- 
bach molecules arises due to the shared dissociation con- 
tinuum |28j . which effectively enhances the weak pho- 
toassociation coupling to be comparable to the Feshbach 
coupling, but it also enhances the Feshbach losses to 
be comparable to the photoassociation losses. In other 
words, the Feshbach molecular state decays vicariously 
through the photoassociation state, and peak enhance- 
ment therefore occurs where the Feshbach detuning is 
large compared to the spontaneous decay rate of the 
photoassociation state, roug hly uj/T s w 5 in Fig. [2} Al- 
though a completely dark state is absent, it so happens 
that the photoassociation state is dark and the Fesh- 
bach state is dim. In weak photoassociation, collisions 
disrupt the dark state and conversion is independent of 
pulse sequence for short pulses, and longer pulses help 
the counterintuitive scheme until the timescale for col- 
lisions is reached. In strong photoassociation, collisions 
play a lesser role. For a counterintuitive pulse sequence, 
the dark photoassociation state then allows efficient con- 
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FIG. 2: (Color online) Feshbach-stimulated Raman pho- 
toassociation of a two-component Bose-Einstein condensate 
for weak (a-d) and strong photoassociation (e-h), for the 
magnetically-tunable (blue dashed, □) and the full quasicon- 
tinuum (red solid, o) models. The dashed (solid) horizontal 
line is the quasicontinuum (tunable-collisions) result for zero 
magnetic field. The dimensionless Feshbach detuning is oj/T s . 



version closer to the Feshbach resonance, and the en- 
hancement peak shifts to the red. For an intuitive pulse 
sequence where the dark state is moot, stronger laser 
coupling requires larger Feshbach detuning to combat vi- 
carious losses, and the peak shifts to the blue. As with 
lasers alone, the peak efficiency decreases with increas- 
ing strength of intra-atomic, molecular, or atom-molecule 
collisions. 

In terms of the effective model, peak improvement 
in the weak case occurs roughly where the resonant 
collisional interaction vanishes Ar = 0, or lu/T s = 
a 2 /(4Ai 2 r s ) sa 9, regardless of pulse sequence. For 
strong photoassociation, collisions again play a lesser 
role, and the counterintuitive scheme is more efficient 



6 



closer to resonance, even closer than in the quasicon- 
tinuum model due to the faux dark state in the effec- 
tive three-level model. For the intuitive pulse sequence, 
the difference is made up-perhaps coincidentally-by the 
peak laser contribution to the resonant collisional in- 
teraction, — S!q/(4<5). For weak or strong photoassocia- 
tion, with tunable collisions nulled, the conversion effi- 
ciency is determined by the non-tunable collisional in- 
teractions. That peak improvement occurs where the 
magnetically-tunable collisional interaction vanishes, and 
that the magnitude of peak improvement is determined 
by non-tunable collisions, is in line with previous work 
on Raman photoassociation of an interacting Bose con- 
densate |I 9, 28: combined with a far-detuned Feshbach 
resonance [28] . So far, so good. 

The role of Feshbach enhancement can be understood 
further by considering the magnetically tunable parame- 
ters in the effective model [Fig[3], where panel (a) corre- 
sponds to the parameters in Fig. [2] The resonant decay 
rate drops relatively quickly to its ambient value, and the 
photoassociation coupling follows shortly thereafter, but 
it is not until the resonant collisional interaction reaches 
a reasonable value (lj/T s < 9) that any improvement 
kicks in. At peak enhancement, the photoassociation 
coupling, weak or strong, is enhanced to FIr/FI = 1.2, 
which amounts to an effective increase in intensity of a 
factor of about f .4. Part of the reason for the lackluster 
improvement in strong photoassociation is due to it al- 
ready being saturated [28 , 3Ql [34] [35] , so any increase in 
coupling is moot. 

The truly unexpected find is that, whereas an in- 
crease in the strength of non-tunable collisions (intra- 
atom, molecule, or atom-molecule) will decrease the 
conversion efficiency [TH] [28J, an increase in the ambi- 
ent value of the tunable inter-atomic collisional interac- 
tion, A12, will actually increase the efficiency of stable 
molecule production. In particular, since the Feshbach 
coupling [If a a -v/A 12 , the Feshbach-detuning location 
of Aj? = 0, i.e., uj z = a 2 /(4Ai 2 ), is independent of Ai 2 . 
At the magnetic-field location of peak enhancement, u) z , 
the resonant contribution to the photoassociation cou- 
pling is an/(2Lj z ) oc A12, and stronger inter-atomic col- 
lisions therefore lead to stronger Feshbach enhancement 
at u> = iij z , as illustrated in Fig. |3jb) for weak photoas- 
sociation and A 12 — > 10A 12 . For the counterintuitive 
scheme at a pulse area of fior = 5 x I0 3 , in the resonant- 
interaction model a factor of two (ten) increase in Ai 2 
enhances conversion from 30% to 37% (60%), and in the 
full quasicontinuum model the former (latter) increase in 
Ai 2 enhances conversion from 30% to 42% (67%). 

These numerical results are supplemented analytically 
as follows. Deriving a two-level system from Eqs. (10) 



in the limit of large Feshbach-shifted detuning {v ^> f^), 
and then deriving a rate equation for the atom losses J2S] 
l3"4"] . we obtain the rate constant for Feshbach-enhanced 
Raman photoassociation for CW lasers 
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FIG. 3: (Color online) Magnetically-tunable parameters vs. 
Feshbach detuning for weak photoassociation, where the solid 
blue line is the photoassociation coupling (J1h/^), the dashed 
red line is the tunable decay rate for the photoassociation 
molecule (Tr/T), and the dot-dashed green line is the tunable 
collisional interaction (At/Ai2). The vertical line denotes 
A_r = 0, the dimensionless Feshbach detuning is defined as 
oj/r a , and the parameters in panel (a) are the same as in 

Fig-i 



where the tunable two-photon coupling is 2\2 = Qrx/ u i 
the tunable two-photon decay rate is 4r 2 = {x/v) 2 Tr, 
and the tunable mean-field shift is approximated to the 
static value cr m f).R = c m f + 2A^ with the non-resonant 
mean field shift a m[ = An + A 22 - (A 13 + A 23 + A 33 ). 
Results shown in Fig. [4] for weak photoassociation with 
v = I0r s are broadly consistent with Fig. [2] In particu- 
lar, the dip below the Feshbach resonance corresponds to 
fin = 0. The peak in atom losses at the Feshbach reso- 
nance corresponds to a peak in two-photon losses, which 
corresponds to the decreased molecule formation near 
Feshbach resonance in Fig. [2] The peak far-above the 
Feshbach resonance arises from Ar = cr m f/2, which cor- 
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FIG. 4: (Color online) Feshbach enhanced two-photon rate 
constant for weak photoassocaition, where the solid (dashed) 
lines corresponds to A12 = 0.04 (0.4). Here the dimensionless 
Feshbach detuning is again defined as ui/Y a , and the dimen- 
sionless rate constant is defined relative to the un-enhanced 
two-photon rate constant for either a non-interacting (Koo) 
or an interacting (Ko) condensate (see text). 



responds to the numerical peak for = 0. We attribute 
the difference to the static mean-field shift approximation 
that neglects transient populations, which also leads to 
a peak location that is not independent of A12. Nev- 
ertheless, for A 12 = 0.04 the enhanced two-photon rate 
constant peaks at roughly the value for a non-interacting 
gas, Koo = X20/T2O) and for A i2 = 0.4 it increases to 
K ~ 5.Koo- Moreover, compared to the a = result 
for an interacting condensate, 4A~ = X^o^o/I^mf^ + 
r| /4], the Feshbach resonance enhances the rate two- 
photon constant by roughly three orders of magnitude 
for A12 = 0.04, which increases to well over four orders 
of magnitude for A 12 = 0.4. Note that 2%2o = ^xl^i 
4r 20 = ( X /S) 2 {T S + 7 pa), and a mf)o = a mf + 2A 12 . 



Before closing, we emphasize that the two-photon de- 
tuning has not been optimized, and that the laser detun- 
ings in general have not been chirped jJT] , both of which 
could lead to further improvements. Also, while compar- 
isons between thermal and condensate systems should 
be taken with a grain of salt, the results for Feshbach 
enhancement in a thermal gas [12] indicate reduced ef- 
ficiency upon averaging over density, and improved effi- 
ciency for narrower Feshbach resonance. Off hand, in as- 
sociation of an interacting condensate the atom-molecule 
coupling oc y/p and the collisional coupling oc p. Colli- 
sions therefore play the biggest role at the center of the 
trap, and inhomogeneity should have little effect on final 
conversion efficiencies in a local density approximation. 
Nevertheless, we look forward to a full investigation- 
including an explicit trapping potential-of condensate 
inhomogeneity for both wide and narrow Feshbach reso- 
nances. 



SUMMARY 



In short, we find that a strong Feshbach resonance can 
substantially improve weak-but not strong-Raman pho- 
toassociation, independent of pulse ordering. The lack of 
Feshbach-enhancement in strong photoassociation is at- 
tributed to an already-saturated transition from atoms to 
molecules, and the independence of pulse ordering to an 
un-optimized pulse length for the counterintuitive pulse 
order. For weak and strong photoassociation, counterin- 
tuitive pulse sequences are indeed more efficient for larger 
pulse areas since the photoassociation molecular state is 
still dark (numerically), even if the Feshbach molecular 
state is dim. 

In the quasicontinuum model, vicarious photoassocia- 
tion losses from the Feshbach state mean that peak en- 
hancement occurs when the Feshbach detuning is large 
compared to the photoassociation linewidth. In the 
resonant-interaction model, peak enhancement occurs 
where the resonant inter-atomic interaction vanishes. 
Also, disagreement between the resonant-interaction and 
quasicontinuum models on final conversion efficiencies 
and the nature of the dark state highlights the impor- 
tance of explicitly including the Feshbach molecular state 
in modeling magnetoassociation. 

Finally, whereas the peak conversion efficiency de- 
creases for stronger intra-atomic, molecular and atom- 
molecule collisions, we find that the peak conversion ef- 
ficiency actually increases for stronger inter-atomic colli- 
sions. Systems with a combination of a strong Feshbach 
resonance and strong inter-atomic collisions will there- 
fore be of greater experimental utility, compared to those 
with a strong Feshbach resonance and weak inter-atomic 
collisions. 
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